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We investigate SU(2) gauge fields topology using new approach, which exploits the well known 
connection between SU(2) gauge theory and quaternionic projective cr-models and allows to formu- 
late the topological charge density entirely in terms of cr-model fields. The method is studied in 
details and for thermalized vacuum configurations is shown to be compatible with overlap-based 
definition. We confirm that the topological charge is distributed in localized four dimensional re- 
gions which, however, are not compatible with instantons. Topological density bulk distribution is 
investigated at different lattice spacings and is shown to possess some universal properties. 
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I. INTRODUCTION 

It is hardly possible to overestimate the significance 
of the topological fluctuations in the vacuum of Yang- 
Mills theories. The discovery of instantons 0, merons Q 
and other topologically non-trivial solutions was the clue 
to understand the major non-perturbative phenomena in 
QCD and to construct rather successful phenomenolog- 
ical low-energy models (see, e.g., Ref. and references 
therein). At the same time the validity of these mod- 
els was in fact questioned from very beginning |4j and 
in particular because in the strong coupling regime they 
become unreliable with no control on the degree of ap- 
proximation made. 

The unique approach to investigate the issue would be 
the numerical lattice simulations. However, until mid 
90's all lattice methods aimed to investigate the topolog- 
ical aspects of the gauge theories were plagued by essen- 
tially the same decease since they could not be applied 
reliably to 'hot' vacuum configurations. Here we mean 
the field-theoretical [f| and geometrical 0, Q definitions 
of the topological charge on the lattice which are known 
to be an art rather then the unambiguous first-principle 
investigation tools (see, e.g., Ref. 8]). The situation was 
ameliorated when the lattice chiral fermions were con- 
structed and a nice alternative topological charge def- 
inition was provided |10| via Atiyah-Singer index theo- 
rem. 

Since then a lot of data had been accumulated and 
we're not in the position to review the current status 
of the problem. Instead we note that strictly speaking 
the usual field-theoretical topological charge and the one 
constructed from fermionic modes are equivalent in the 
continuum limit and for smooth fields only. In the con- 
text of regularized quantum field theory they are a priori 
distinct and are sensitive quite differently to the details 
of regularization. It is believed that the continuum limit 
of both definitions is the same, but at finite cutoff the 
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overlap-based construction is much more advantageous 
although it is not ideal either, in particular, because of 
its complexity. 

The purpose of this paper is to develop an alternative 
approach to the topology of D — 4 pure SU(2) lattice 
gauge theory, which ideologically is very close to what 
had been said above. Namely, we propose to consider the 
definition of the topological charge which is known to be 
equivalent to the standard one in the continuum limit. 
However, at finite cutoff it becomes an independent con- 
struction and for this reason is worth to be investigated. 
Our approach could be best illustrated in the context 
of D = 2 CP 1 (or 0(3)) cr-model where the topological 
charge 

Q CP i = — I d 2 xe^ u n ■ [djt x d v n] (1) 

87T J 

could be identically rewritten in terms of U(l) gauge po- 
tentials Aft — —iz*dfj,z 

Qcpi =^Jd 2 xe^F UVl (2) 

where is the usual Abelian field-strength and com- 
plex z is related to n by standard stereographic projec- 
tion. The idea is to follow the above equations in the 
opposite order and formulate the topological charge of 
D = 2 U(l) gauge fields as the corresponding one in 
CP™ cr-model, where we have indicated that in fact the 
cr-model rank n is a free parameter. The above reasoning 
generalizes almost trivially to SU(2) case. The difference 
is that we have to consider the field of real quaternions 
H instead of the complex numbers and the corresponding 
D = 4 cr-model target space is the quaternionic projec- 
tive space HP™. Note that the deep connection between 
HP™ cr-models and SU(2) Yang-Mills theory is in fact 
well known (see, e.g., Refs. Il2l and references 
therein) . In particular, the geometry of gauge fields could 
best be analyzed in the HP" context. In fact, all known 
instantonic solutions of Yang-Mills theory could be in- 
duced from the topological configurations of suitable HP™ 
cr-model. 
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Thus it is natural to expect that the construction of 
'nearest' to the given gauge background HP" fields cap- 
tures accurately the gauge fields topology leaving aside 
their non-topological properties. Clearly, the notion of 
'nearest' is crucial and it is described in details below. 
Here we note that it inevitably introduces a sort of ambi- 
guity in our approach the significance of which is hardly 
possible to analyze a priori for quantum vacuum con- 
figurations. Although for classical gauge fields the cor- 
responding HP" cr-model is known to be unique, the 
mathematical rigour is lost in case of 'hot' vacuum fields 
and the corresponding systematic errors are to be in- 
vestigated separately, presumably using numerical meth- 
ods. This concerns, in particular, the topological charge 
defined via HP" cr-model embedding and therefore its 
properties near the continuum limit are to be carefully 
studied. These issues are considered only preliminary in 
present publication and therefore the values of physical 
quantities quoted below are to be taken with care. Note, 
however, that our results do indicate that the ambigu- 
ity of HP" cr-model construction for equilibrium vacuum 
configurations is likely to be irrelevant while other sys- 
tematic errors could easily be reduced algorithmically. 

The paper is organized as follows. In section [H] we 
give basic theoretical background, introduce our nota- 
tions and describe the proposed approach in the contin- 
uum terms. Section IIIII is devoted to the detailed de- 
scription and investigation of our method on the lattice. 
In particular, we consider the dependence of the topo- 
logical charge and its density definitions upon the HP" 
cr-models rank and other parameters involved. Then the 
algorithm is tested on semiclassical fields and compared 
with overlap-based topological charge on 'hot' SU(2) con- 
figurations, where we found a quite remarkable agree- 
ment. 

Section llVl is devoted exclusively to the investigation 
of SU(2) gauge fields topology as it is seen by our method 
and a few closely related issues. We show in section HV Al 
that the dynamics of embedded HP 1 cr-model closely re- 
flects the dynamics of the original gauge fields, in par- 
ticular, the mass gap of the cr-model is given by SU(2) 
string tension. Next in section llV Bl the topological sus- 
ceptibility is considered and we find that it scales rather 
precisely being in agreement with the existing literature. 
Then in section ITvTI we investigate in details the lo- 
cal structure of the topological charge at the particu- 
lar value 0.1193(9) fm of the lattice spacing. We con- 
firm that the topological charge bulk distribution has 
rather pecu liar lumpy structure discovered in previous 
studies [Til ITEl ITil 1 1 7j which, however, has nothing to do 
with instantons. In particular, the lumps are distributed 
according to N(Vc) ~ V c 3 ^ 2 , where Vc is the lump 4- 
volume, and are dominated by UV small lumps. Nev- 
ertheless, the topological susceptibility turns out to be 
lumps saturated. It is amusing that setting the param- 
eters involved in the lumps definition to unphysical (as 
we believe) values we are able to qualitatively reproduce 
the recently discovered 0] global topological structures 



( section llV C 2|) . Namely, the lumps organize themselves 
into a few (typically two) percolating structures per con- 
figuration plus the divergent amount of UV small lumps 
consisting mainly of just one lattice point. However, for 
reasons which we discuss in details it is unclear for us 
what is the physical significance of this result. Then in 
section llV C 3l the topological density correlation function 
is considered and we argue that it should not be neces- 
sary negative within our approach. Instead it reflects the 
lumpy structure mentioned above and the correspond- 
ing correlation length is determined by the characteristic 
lump size, which turns out to be of order half the lat- 
tice spacing. Finally, in section llV C 41 we confront the 
data obtained at two different f3 values and argue that 
the parameters involved in the lumps definition are likely 
to be given in physical units. Moreover, the lumps vol- 
ume distribution N(Vc) ~ U^T 3 ^ 2 seems to be spacing 
independent. As far as the characteristic lump size is 
concerned, we are not confident in its scaling properties, 
more data is needed to quantify the issue. 



II. HP™ d-MODELS AND THEIR EMBEDDING 
INTO SU(2) YANG-MILLS THEORY 

The purpose of this section is to outline in the con- 
tinuum terms the approach we proposed to investigate 
SU(2) gauge fields topology. In particular, we briefly re- 
mind the construction of quaternionic projective spaces 
HP" and corresponding four dimensional cr-models. The 
material of this section is by no means new, for review 
and more details see, e.g., Refs. [H[ll [II. 



A. Quaternionic Projective Spaces and HP™ 
cr-models 



HP" is an example of quaternionic Grassmann man- 
ifold and is a compact symmetric space of quaternionic 
dimensionality n on which the group Sp(n + 1) acts tran- 
sitively. It can be viewed as the factor space 



HP r 



Sp{n + \) 



(3) 



Sp{n) x Sp(l) 

and therefore for arbitrary n we have 

^ 4 (HP") - 7r 3 (Sp(l)) = n 3 (SU(2)) = Z. (4) 

Another useful representation of HP" is 

HP" = S 4n+3 /S 3 , (5) 

hence HP" is the fibering of S 4n+3 over quaternionic lines 
passing through origin. In the particular case n = 1, 
which will be the most important below, we have 



HP 1 = S 7 /S 3 



(6) 
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(second Hopf fibering). The simplest explicit 
parametrization of HP" is as follows. Consider normal- 
ized quaternionic vectors 

n 

\q) = [qo,--,q n } T , (?l9>=^i?i = ieH, (7) 

where [...] T denotes transposition, = qf e a , a = 0, 3 
are real quaternions (homogeneous coordinates) with e a 
being the quaternionic units 



eo = 1 , 



-Sij+eijkek, i,j,k = 1,2,3 (8) 



and quaternion conjugation is defined by 

q~i=q?e a , e = e , ej = -ei, i= 1,2,3. (9) 

Therefore the states | q ) describe the (4n+3)-dimensional 
sphere, \q) G S 4 "+ 3 . According to Eq. © the HP™ 
space is the set of equivalence classes of vectors | q ) with 
respect to the right multiplication by unit quaternions 
(right action of SU(2) gauge group) 



Q 



\q)v, 



1, v G H. 



(10) 



As usual it is convenient to introduce quaternionic (n 
1) x (n + 1) projection matrices 



P = \q)(q\ (P lk =q t q k ), 



P 



P = pt, (11) 



which are invariant under IjlUfl and therefore parametrize 
HP™. Alternatively, one could consider the matrix 



N = 2P - 1 



N = N*>, N 2 



1 , 



(12) 



which generalizes to HP™ case the familiar stereographic 
projection of S 4 = HP 1 . Indeed, for n = 1 the two- 
component vector (J7J is uniquely defined by inhomoge- 
neous coordinate u> — qiq^ 1 G H and in terms of unit five 



dimensional vector n 
2lu 



A 



n 



1 



a = 0, 3 , 



(13) 



we have N = Xm=o l AnA ? where j A are the five Eu- 
clidean Dirac matrices {7™,7 5 }. 

As far as the HP™ cr-models are concerned it is clear 
from Eq. Q that the model is non-trivial provided that 
the base space is taken to be the 4-sphere S 4 . The in- 
dex of the mapping S 4 — ► HP™ is given in terms of the 
projectors fTT|l by 

Q = ^ J d 4 x Sctr (Pd^Pd u Pd x Pd p P) , (14) 

where tr means the trace over quaternionic indices and 
scalar part is defined by Scq = (q + q)/2. The action of 
the model is given also in terms of projectors ljll|l . but 
its explicit form will not be needed in what follows. Note 
that Eq. (|14|) simplifies greatly in n = 1 case. In terms of 



five dimensional unit vector (|13f) the topological charge 
is 

Q = T^yj J d 4 x e^ Xp £ abcde • (15) 

• n A d„n B d u n c d x n D d p n E 

and geometrically is the sum of the oriented infinitesimal 
volumes in the image of the mapping n A (x) : S 4 — > S 4 . 

An alternative way to describe the HP™ cr-models is 
to consider auxiliary SU(2) gauge fields transforming 
as usual 
derivative 



vA^v — vd^v under (|10f> . The covariant 



D\q)=d\q) + \q)A 



(16) 



transforms homogeneously D\q) — > D\q)v and both 
the topological charge and the action could be written 
equally in terms of D \ q). The potentials A^ are not dy- 
namical degrees of freedom and could be eliminated in 
favor of | q ) 



A„ 



[q\d^\q) 



,qSu,qi 



i=0 



(17) 



However, it is crucial that the topological charge H14fl . 
(|15ll is expressible solely in terms of A^ 

Q = sk^J d4x £ ^ Xp Tr F ^ Fxp ' (18) 

Ffiv — dfiA v — d y A^ + [Afj,, A u ] , 

being essentially equivalent to the familiar topological 
charge of the gauge fields {T7J|. A particular illustration 
is provided by hedge-hog configuration of n A fields of 
HP 1 cr-model. One can show that indeed Eq. I|17(l leads 
in this case to the canonical instanton solution of SU(2) 
Yang-Mills theory. 

B. Relation to SU(2) Yang-Mills Theory 

The difference between the HP™ cr-models and the 
usual SU(2) Yang-Mills theory lies in the fact the gauge 
potentials l|17|) are not independent but composite fields 
constructed in terms of elementary scalars. However, this 
difference becomes less significant due to the theorem by 
Narasimhan and Ramanan ^| (see also [l3|). Namely, 
it was shown that any sourceless SU(2) gauge fields con- 
sidered on a compact manifold can be induced from the 
canonical connection of suitable Stiefel bundle (it goes 
without saying that the theorem applies to smooth fields 
only). In particular, that means that any classical SU(2) 
gauge potentials on S 4 could be realized it terms of HP™ 
CT-ficlds as in Eq. (|17fl for sufficiently large but finite n. 
It should be noted that in fact the representation (I17[) is 
well known and was used, in particular, in Ref. |20j to 
construct the general instanton solutions in Yang-Mills 
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theory. Moreover, one could formulate [21j the gauge 
theories entirely in terms of projectors 

Notice also a particular logic underlying the considera- 
tions above. Namely, one starts from the HP™ cr-fields the 
topology of which could be analyzed either directly with 
Eq. (|14fl or indirectly with Eq. (|18fl when the gauge fields 
(|17J) are introduced. In the continuum limit once the dif- 
ferentiability is assumed these two approaches are essen- 
tially the same. However, it is crucial that in quantum 
field theory regularized with finite UV cutoff [21j these 
two approaches need not be identical anymore. More- 
over at finite cutoff various methods are expected to be 
sensitive quite differently to the details of regularization 
and, in particular, to UV-scale fluctuations and lattice 
artifacts. The well known example of this kind is pro- 
vided by the topological charge l|18l) in Yang-Mills the- 
ory for which the field-theoretical [5| , geom etrical 0, QJ 
and Atiyah-Singer index theorem based [jlJ,|22j construc- 
tions are quite distinct at finite a (for review and further 
references see, e.g., Ref. |23|). 

Since our ultimate goal is to investigate the topology 
of Yang-Mills fields at finite UV cutoff, we could try to 
reverse the above logic and attempt to use Eqs. lfT4")l . i(T5|) 
in the context of gauge theories. Therefore generically 
our approach looks as follows. 

i) For given SU(2) gauge potentials and for given rank 
n one finds the corresponding closest HP" cr-model fields 
for which the approximation (fTTj) is as best as possible. 
The quality of approximation a priori depends upon the 
rank and deserve a separate investigation. 

ii) Then the topology of the gauge fields could be 
analyzed in terms of HP" cr-model. In particular, the 
topological charge density is given by oriented four- 
dimensional volume of infinitesimal tetrahedron in the 
image of the map P : S 4 — > HP" provided by projectors 
CD- 

It is clear that the actual implementation of the above 
program is challenging, there are a wealth of technical 
issues to be discussed below. However, right now let us 
mention the following. 

a) Comparing the number of degrees of freedom one 
concludes that the approximation (|17|l could not be ex- 
act for small rank n. At the same time it is well known 
that the classical instanton solution in Yang-Mills theory 
is exactly reproduced by HP" cr-model already for n = 1 . 
This means, in particular, that for gauge fields represent- 
ing classical instanton plus perturbative fluctuation the 
above construction for n — 1 is likely to reproduce cor- 
rectly the gauge fields topology while being almost blind 
to the perturbative noise. Naively, one could expect that 
the approximation (|17|l is to be more stringent as rank is 
increased. Then it follows that the rank of the HP" cr- 
model could play the role of ultraviolet filter: the larger 
the rank n the better is the approximation l|17|) and there- 
fore the above approach becomes more sensitive to the 
ultraviolet fluctuations. 

b) In the continuum limit the topological charge den- 
sity is given by the oriented volume of infinitesimal 



4-dimensional tetrahedron embedded into HP" space. 
However, at any small but finite lattice spacing the tetra- 
hedron need not be infinitesimal and its vertices could be 
far away from each other (in the standard HP™ metric). 
The commonly accepted prescription in this case (see, 
e.g., Refs. 0,0 0i 113) is to connect vertices pairwise 
by shortest geodesies and calculate the 4- volume of corre- 
sponding 'geodesic' tetrahedron. Unfortunately, it seems 
that this procedure is impossible to implements practi- 
cally for n > 1. Indeed, to the best of our knowledge 
even in the simplest case n = 1 the analytical expres- 
sion for the volume of 4-dimensional spherical tetrahe- 
dron embedded into S 4 space does not exist. Moreover, 
for n > 1 we don't see any practical way to calculate the 
volume even numerically. However, for n = 1 the simple 
geometry of the space HP 1 = S 4 admits a reasonable nu- 
merical solution (see below for details). Therefore in the 
investigation of the gauge fields topology we are forced 
to consider the HP 1 cr-model only. However, in view of 
the item a) above this should not be considered as fatal 
restriction of our method (see also the next section). 



III. LATTICE IMPLEMENTATION 

In this section we describe in details the implementa- 
tion of the above approach on the lattice. Section llll Al is 
devoted to the problem of cr-model embedding for given 
SU(2) lattice gauge fields background. In particular, 
we discuss the dependence of approximation (|17|l upon 
the rank n and investigate the Gribov copies issue. In 
section IIII Bl we describe the topological charge density 
calculation method which uses solely the HP 1 cr-model 
fields. 



A. HP" cr-models Embedding on the Lattice 

The construction of HP" cr-model starts from assigning 
quaternionic vectors \q x ), Eq. Q) , to each lattice site and 
then the actual problem is how Eq. i|17|) translates to the 
lattice. The l.h.s. of Eq. I|17fl should evidently be given 
by link matrices U x .^. As far as the r.h.s. is concerned 
its lattice counterpart is in principle ambiguous. In fact, 
the correct lattice replacement of Eq. (fT7|) could be found 
by analogy with corresponding case of CP" cr-models [24| 
or deduced from considerations of quaternionic quantum 
mechanics 261 



[ Qx | Qx+n > 



(19) 



\(q x | q x +n)\ ' 

Since generically it is not possible to satisfy Eq. i|19|) iden- 
tically the best what we can do is to minimize the func- 
tional 



X 

5F x = ±£(l-Sc[£/t 



(20) 



4i) 



5 



x- x- 



HP 1 
HP 8 
HP 15 



7000 - 
6000 - 
5000 - 
4000 - 
3000 - 
2000 - 
1000 - 



X 

X x 



o 



+ x- ^X-X X ^ 



+ x 
+ 

X + 

X- 



Cooled conf. 1 + 
Cooled conf. 2 
Cooled conf. 3 X 



FIG. 1 

12 4 , /3 



: Typical distribution of SF X , Eq. 12U1 , on thermalized 
= 2.35 lattices for various cr-model ranks. 



FIG. 2: Distribution of 5F X , Eq. $ZU}, for HP 1 cr-model on 
three highly cooled 12 4 configurations containing respectively 
three instantons, one anti-instanton and no instantons at all. 



over all {| q x )} configurations (here V denotes the lattice 
volume). Note that the minimization tasks of this kind 
are expected to be plagued by Gribov copies problem to 
be discussed below. We implemented the usual local over- 
relaxed minimization of the functional l|20() . The stop- 
ping criterion is that the difference of F values between 
two consecutive sweeps should be smaller then 10~ 6 . 

The first question to be addressed is the quality of the 
approximation Ijl9(l , the natural measure of which is given 
by (SF X ). We measured (SF X ) on about 20 thermalized 
statistically independent SU(2) gauge configurations on 
12 4 lattice with (3 = 2.35 (the rank of the cr-model was 
n = 1). The averaged value turns out to be (SF X ) w 
0.12 (see also below) which in fact indicates rather good 
approximation quality. What is more important is that 
(SF X ) appears to be slightly decreasing function of (3. For 
instance, on 16 4 lattice at (3 = 2.40 we found (SF X ) ss 
0.10. 

Next we have to investigate how the quality of ap- 
proximation depends upon the cr-model rank n. We cal- 
culated the distribution of 6F X for n — 1,8,15 on sev- 
eral fixed thermalized SU(2) configurations (12 4 lattice 
at f3 = 2.35). The corresponding typical histograms are 
presented on Fig. and clearly show that the approx- 
imation quality does not depend at all upon the rank. 
Thus the naive expectation that rank could serve as an 
ultraviolet filter turns out to be wrong. Instead we found 
that the approximation quality is 10-20% for 'hot' SU(2) 
gauge fields irrespectively of the concrete HP" cr-model 
considered. The plausible explanation of this might be as 
follows: from the continuum considerations we know that 
the equality 117|l is valid for smooth (e.g., classical) SU(2) 
fields only. It is commonly believed that the configura- 
tions at (3 — 2.35 are not smooth. For this reason the 
approximation (|19fl could not be exact for any rank and 
its quality is mainly given by the gauge fields roughness 
and indeed rank independent. Note that if this reason- 
ing is valid then the approximation quality l|20|) provides 
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FIG. 3: SF X distribution on 12 4 lattice for random input gauge 
fields. 



a quantitative measure of gauge fields smoothness. The 
conclusion is that it is unnecessary to consider HP™ a- 
models with n > 1 and from now on we confine ourselves 
to the case n = 1. 

In view of the above discussion it is instructive to con- 
sider semiclassical and completely random inp ut g auge 
fields. Using the cooling procedure (see, e.g., [23 and 
references therein) we generated three highly cooled 12 4 
configurations containing respectively three instantons, 
one anti-instanton and no instantons at all. The corre- 
sponding distribution of SF X is shown on Fig. It is 
worth to mention that the approximation quality l|19(l 
turns out be excellent, the gauge fields are reproduced 
with accuracy < 1% in the whole lattice volume. Note 
that for non-trivial classical background the distribution 
of 8F X is slightly broader then that for trivial case still 
reproducing the input gauge fields almost exactly. As 
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(with mean value subtracted) calculated on 12 4 (3 = 2.35 
thermalized lattices and averaged over 20 configurations, see 
text. 



far as the random gauge fields are concerned, the re- 
sults obtained in this case (Fig. |3J) support the above 
argumentation. Namely, the distribution of 8F X values is 
much broader being significantly large even at 8F X « 0.5. 
The corresponding quality of approximation could be es- 
timated as {8 F x ) ~ 0.3 and indeed is much worse then 
that in previous cases. 

The last problem to be considered in this section is the 
problem of Gribov copies usually associated with mini- 
mization of functionals like (|20|l . The problem lies in 
the fact that there might be many local minima of 120(1 
which are relatively far apart in the configuration space 
of | q x ) fields. To investigate this issue we performed the 
following calculation. For a particular thermalized gauge 
configuration (12 4 lattice at (3 = 2.35) the minimization 
of H20fl was repeated 50 times starting every time with 
some random \q x ). From the deviations of found min- 
imal values from their mean the distribution histogram 
was constructed and averaged over 20 input gauge config- 
urations. In order to quantitatively address the Gribov 
copies problem one should compare the width of F val- 
ues distribution obtained this way and the precision of 
the functional l|20|l calculation which in our case could 
be estimated as ~ 10~ 6 . The distribution of the minimal 
values of F is shown on Fig. 01 and its width turns out to 
be of order i=s 10 -5 . The conclusion is that the Gribov 
copies problem is likely to be inessential in finding the 
HP 1 er-model for given SU(2) gauge background. How- 
ever, for safety reasons we always tried 10 random initial 
configurations | q x ) in our minimization procedure and 
then selected the best minimum found. 



B. Topological Charge Density 

As we already mentioned the topological charge den- 
sity in terms of the HP 1 cr-model fields is given by the 
oriented 4-volume of spherical tetrahedron T embedded 
into S . The necessity to consider 4-dimensional sim- 
plices could be deduced already from the continuum ex- 
pression (|15|) which depends upon unit vectors n^k at five 
infinitesimally close neighboring points marked by index 
i = 0, ...,4. On the lattice the vertices of the tetrahedron 
T = {n^} G S 4 become finitely separated and should be 
connected pairwise by shortest geodesies thus forming a 
spherical tetrahedron embedded into S 4 . 

To the best of our knowledge the volume of spher- 
ical four dimensional tetrahedron is unknown analyti- 
cally contrary to the case of three dimensions (see, e.g., 
Refs. j2^| and references therein). Therefore, the only 
way to proceed is to invent some reasonable numerical 
algorithm to evaluate the 4-volume in question. Note 
that this is so because we are specifically interested in 
the topological charge density. As far as only the global 
topological charge is of concern, one could proceed in 
usual way utilizing the fact that the topological density is 
almost the total derivative (see, e.g., Ref. 25] ). Then the 
evaluation of the global charge reduces to the calculation 
of the volumes of three-dimensional spherical simpliccs 
which could be done analytically. 

The most straightforward way to estimate the 4- 
volume of spherical simplex is to use the Monte Carlo 
technique, the essence of which is to pick up a random 
point in S and to decide whether it is inside or outside 
the spherical tetrahedron. Due to the simplicity of S 4 
geometry the last question reduces to the relatively fast 
calculation 5x5 determinants analogously to the case 
of S 3 sphere considered in 0. However, proceeding this 
way we definitely loose the integer valuedness of the topo- 
logical charge since volume of every simplex (topological 
charge density) is calculated with finite accuracy. Evi- 
dently the accuracy should be small enough to recover 
the topological charge and simultaneously large enough 
to be practically reasonable. In our test runs (see below) 
we found that the optimal relative accuracy of the spheri- 
cal tetrahedron volume estimation should be smaller then 
2%. 

Now we can summarize our algorithm of topological 
charge density calculation. 

1. On input there are five unit 5-dimensional vectors 
n (i)' i = 0, ...,4 sitting at five vertices of simplex T of the 
physical space triangulation. The vectors are con- 
structed via Eq. (|13|) and provide the mapping of T into 
the spherical tetrahedron T 6 S 4 , the oriented volume 
V(T) of which is to be calculated. Generically the set 

is not degenerated det^ [nfy] ^ 0. Note that for 

|detiA [nA]\ < 10~ 10 we equate V(T) to zero. 

2. In order to speed up the Monte Carlo integration 
we calculate the normalized mean vector N oc 

and find the 5-dimensional cone C around N containing 
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FIG. 5: Distribution of the topological charge Q float, 
Eq. 1221 . for single 5-dimensional simplex with random HP 1 
cr-model fields at its vertices. 

all the n {i) . Since V(T) < V(C n S 4 ) the tetrahedron 
volume is equated to zero if V(C fl S 4 ) < 10~ 8 . 

3. The Monte Carlo estimation of V(T) begins with 
picking up 10 3 random points uniformly distributed in 
CDS 4 . If none of them falls into the interior of T the 
tetrahedron volume is equated to zero. Otherwise we 
continue to generate uniformly in C n S 4 random points 
until the volume V(T) is known with accuracy 2%. Fi- 
nally, the topological charge density in the simplex T of 
physical space triangulation is given by 

q(T) = JL sign (det lA [n&]) • V(T) , (21) 

where the normalization factor 87r 2 /3 is the volume of 
S 4 . The corresponding topological charge is 

Q/Ioat = $>CO , (22) 
T 

where we have explicitly indicated that it is not integer 
valued. We remark that in order exploit this procedure 
on usual hypercubical lattices the hypercubes should be 
sliced into simplices. We did this according to the pre- 
scription of Ref. j^jj . 

The disadvantage of the algorithm is that the topo- 
logical charge is not integer valued. However, for small 
enough accuracy of calculation at each simplex the prob- 
lem is likely to be inessential. In order to check this we 
performed the following tests. First, we generated 500 
sets of six random 5-dimensional unit vectors which are 
geometrically assigned to the vertices of 5-dimensional 
tetrahedron, the boundary of which consists of six 4- 
dimensional simplices and is S 4 topologically. Each gen- 
erated set is the mapping S 4 — > S 4 and we applied our 
algorithm to calculate the distribution of Q float- The re- 
sulting histogram is presented on Fig.[S]and clearly shows 
that for this test the chosen accuracy is more then enough 
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FIG. 6: Distribution of Q float- Q, Eqs. calculated 
on 55 configurations thermalized at /3 = 2.40 on 16 4 lattice. 

since the topological charge is sharply peaked around 
integer values. However, this simple test could not be 
completely convincing since in Eq. 1221) the uncertainties 
coming from each simplex T are accumulating. Thus 
the deviation of the topological charge 1)221) from near- 
est integer increases with lattice volume V and for given 
accuracy per each simplex the distribution of the topo- 
logical charge flattens in the limit V — > 00. Therefore 
we must ensure that for the lattices we're going to con- 
sider and for the chosen accuracy per each simplex the 
distribution of Qfi oa t is indeed peaked around integers. 
Below we discuss the physically relevant results obtained 
on 16 4 lattices at j3 = 2.4, 2.475. Since the distribution 
of Q float is technical rather then physical issue it is pre- 
sented in this section. The distribution of non-integer 
part of the topological charge 122[) calculated on 16 4 lat- 
tice at (3 = 2.40 with 55 configurations is shown on Fig. [H] 
As expected the histogram is much broader then that on 
Fig. El but nevertheless is still reasonably sharp with dis- 
tribution width being w 0.2. We conclude therefore that 
for the chosen parameters and lattice volumes we could 
safely identify the topological charge as 

Q= [Qfioat], (23) 

where [x] denotes the nearest to x integer number. We 
stress that the identification (|23|l is not universal and is 
valid only for particular range of parameters and lattice 
volumes. From now on the term "topological charge" will 
refer exclusively to the integer valued quantity Q. 

Unfortunately, the presented algorithm turns out to 
be very time consuming. As a matter of fact there is 
no hope to use it on single processor machines. Indeed, 
the considered triangulation of one 4-dimensional hyper- 
cube consists of 16 simplices and their total number be- 
comes prohibitively large for any reasonable lattice vol- 
umes. However, it is very easy to implement the method 
on computer clusters since calculation of the topological 
charge density is done independently for each simplex. 
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FIG. 7: Time slice of the action and the topological charge 
densities 1211 on cooled 12 4 configuration containing three 
instantons. 



FIG. 9: Correlation function 12611 of the overlap- and HP 1 
(j-model based definitions of the topological charge density, 
sec text. 
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FIG. 8: Time slice of the action and the topological charge 
densities 12 1 1 on cooled 12 4 configuration containing one anti- 
instanton. 



1. Testing the Algorithm 

The approximate integer valuedness of the topological 
charge is not the only test which our algorithm should 
pass. At least we should convince ourselves that the 
method works in the case of (quasi)classical configura- 
tions and compare it with known topological charge con- 
structions. 

As far as the quasi-classical fields are concerned we 
applied our method to the same set of highly cooled 12 4 
configurations generated in section IlII Al We have seen 
that the approximation quality of SU(2) gauge matrices 
by the corresponding HP 1 a- model fields was almost ex- 
cellent (about 1%). Now we compare the local gauge 
action and the density of the topological charge on these 
configurations. It turns out that they agree just nicely 



as shown on Fig. \7\ [SJ where for a particular time slice 
t = 2 the amount of red (green) points per 3-dimensional 
lattice cube is proportional to the action (absolute value 
of the topological charge) density. On the first configu- 
ration (Fig. 0) our algorithm identified three instantons 
with total topological charge being Q float — 2.91. On the 
second configuration (Fig. [HJ we found one anti-instanton 
and the topological charge is Qfi oa t = —0.98. Note that 
the regions where the topological charge density is sig- 
nificant are rather large, their volumes are about 3 4 in 
lattice units. 

To summarize, we are confident that our method works 
as expected on quasi-classical configurations. Next we 
would like to compare it with becoming standard nowa- 
days overlap-based topological charge definition 0] . 



2. Comparison with Overlap-Based Topological Charge 

The validity of the topological charge construction we 
discussed so far could only be convincing provided that 
we confront it with other known topological charge defini- 
tions. As far as field-theoretical p| and geometrical 0,0 
constructions are concerned, they could be applied in fact 
only on cooled configurations and we refuse the corre- 
sponding comparison for this reason. On the other hand, 
the definition |To| based on the overlap Dirac operator [t| 



q (x) = -tr 7 5 \1--D x , x j (24) 

is free from the above disadvantages. For reasons to be 
discussed in section FtV C 31 we will not consider the full 
trace in Eq. I|24|) . but restrict it to N lowest eigenmodes 
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a,fm 


L 3 


L t 




Nconf 


2.40 


0.1193(9) 


16 


16 


13.3(4) 


55 


2.475 


0.0913(6) 


16 


16 


4.6(1) 


18 



TABLE I: Simulation parameters. 



= AiV (i) 



^(*) = -£(i-y)c {4) (*)> 

i<JV 

c«(x) = ^ )t 7 5 ^ ) - 



(25) 



The quantity q^ is known as the effective topological 
density [30| . 

In order to make a quantitative comparison of the def- 
initions (|21|l and lll'ol) we measured the correlation func- 
tion |3lJ 



(26) 



on 18 statistically independent 16 4 configurations ther- 
malized at (3 = 2.475. The results corresponding to 
N = 10 and N = 20 lowest eigenmodes taken into ac- 
count are presented on Fig. [§] It is apparent that both 
definitions are indeed compatible with each other and de- 
termine the topological density consistently. Moreover, 
the graph shown on Fig. El is very similar to the corre- 
sponding one given in Ref. [31|, where the overlap-based 
and field-theoretical constructions were compared. 

To summarize, we found clear evidences that our con- 
struction of the topological charge density is locally cor- 
related with overlap-based definition. This provides a 
stringent test of our approach and allows us move directly 
to the investigation of the SU(2) gauge fields topology. 



IV. TOPOLOGY OF SU(2) GAUGE FIELDS 

In this section we discuss the measurements performed 
with HP 1 cr-model construction of the topological charge 
in pure SU(2) lattice gauge theory. In section ITV Al the 
dynamics of embedded HP 1 cr-model is briefly consid- 
ered. Then in section IIVBI we present our results con- 
cerning global topological charge distribution and topo- 
logical susceptibility. Section llV CI is devoted to the de- 
tailed investigation of the topological charge density bulk 
distribution and the corresponding correlation function. 
Our calculations were performed on two sets (Table P) of 
statistically independent SU(2) gauge configurations gen- 
erated with standard Wilson action. The lattice spa cing 
values quoted in the Table [I] are taken from Ref. [32] 
and fixed by the physical value of SU(2) string tension 
\J~o = 440 MeV. As we noted already the method we em- 
ployed for topological charge density calculation is very 
time consuming and this explains the relatively small 
number of configurations we were able to analyze. 
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FIG. 10: The correlation function (12711 in the embedded HP 1 
a- model. Lines correspond to the best fitting curves 1281 . 



A. Dynamics of the Embedded HP 1 a-model 

Given that in our approach the calculation of the topo- 
logical properties of SU(2) gauge fields requires the con- 
struction of HP 1 cr-model it is natural to consider the 
dynamics of the embedded scalar fields. The basic corre- 
lation function in the HP™ cr-models is given by 



G{x) = (Sctr[P(0)P(x)]> 



(27) 



where Pik are the quaternionic valued projector matrices 
ifTTJl. Note that we didn't subtracted the disconnected 
part of G(x). The physical meaning of the correlation 
function l|27l) could well be illustrated for HP 1 cr-model 
where it reduces up to irrelevant constant contributions 
to G(x) = {Y^,A n o n x)- On general grounds one expects 
that the embedded HP" cr-models could be in two dif- 
ferent phases distinguished by G(x). Namely, the corre- 
lator (|27ll falls off exponentially at large \x\ in massive 
(disordered) phase and only with some inverse power of 
\x\ in the massless (ordered) phase. Taking into account 
the quality of the approximation (|19fl discussed in sec- 
tion ^^21 one could argue that the confinement (decon- 
finement) phase of gluodynamics might be reflected by 
the disordered (ordered) phase in the corresponding cr- 
models. At present it is too early to discuss this conjec- 
ture quantitatively, but our preliminary studies indicate 
that it is indeed reasonable. The data we have right now 
refer to the confinement phase exclusively where we ex- 
pect the exponential behavior of the correlator (|27|l 



G(\o 



Co 



i\x\ 



Ci 



(28) 



The results of our measurements of the correlation 
function l|27|) are presented on Fig. 1101 (left). To inves- 
tigate the large distance behavior of G{x) we fitted the 
data to Eq. I|28|l in various ranges of R = \x\. The fits 
turn out to be rather stable, the best fitting curves are 
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FIG. 11: Topological charge distribution calculated on j3 — 
2.40 configurations. 



shown on Fig. El by dashed lines. The results of vari- 
ous fits strongly suggest that the parameter m entering 
Eq. i(2"5|) is nothing but the SU(2) string tension 

m = asfa, (29) 

being equal to m,g = 2.4o = 0.28(2) and 771/3=2.475 = 0.19(1) 
at the (3 values considered. In fact, Eq. indirectly 
supports the above made conjecture on the phase dia- 
gram of the embedded HP 1 a- model. 



B. Topological Charge 

We should mention that the limited set of configura- 
tion we have in our disposal prevents us from studying 
the topological charge distribution and related quantities 
with sufficient accuracy. Nevertheless our results allow 
to draw some preliminary conclusions about the method 
we presented. There are two basic quantities to be dis- 
cussed here, the topological charge l|23|) distribution and 
the topological susceptibility \ = (Q 2 )/V- 

The distribution of the topological charge for (3 = 2.40 
is shown on Fig. EI Although the width of the histogram 
is roughly what is expected at this (3 value, the shape of 
the distribution seems to indicate the lack of statistics. 
In particular, it is hardly possible to fit it with Gaussian 
profile. Moreover, the noticeable asymmetry in odd/even 
Q values should also be mentioned, which probably hints 
on some subtle issues in our topological charge defini- 
tion. Anyhow, the evident lack of statistics prohibits 
us to make solid conclusions at the moment. For the 
same reason we do not present the analogous histogram 
at p = 2.475 since the statistics in this case is even worse. 

As far as the (Q 2 ) and the topological susceptibility 




FIG. 12: Visualization of the topological charge density distri- 
bution in a particular time slice of 16 4 , f3 = 2.40 thermalized 
configuration. 



are concerned, their values turn out to be 



p 


(Q 2 ) 


X 


2.40 


10.8 ± 1.8 


(190(10) MeV) 4 


2.475 


3.9 ±1.0 


(192(12) MeV) 4 



which indicate a perfect scaling and are fairly consistent 
with data of Ref. H2. 

To summarize, at present our calculation of the global 
topological charge is plagued by the lack of statistics and 
should be improved to quantify our method. Neverthe- 
less, qualitatively our results are in agreement with the 
literature and don't show any pathology in the approach 
we have developed. 

C. Topological Charge Density 

The topological charge density q{x) is defined by 
Eq. H21J) and provides quantitative characterization the 
topological charge bulk distribution. As a first qualita- 
tive illustration we present visualization of q(x) in par- 
ticular time-slice of four dimensional lattice. It turns 
out that the corresponding pictures at various slices look 
very similar. The typical distribution of q{x) is shown on 
Fig. U3 (data is taken on 16 4 , p = 2.40 lattices). Here 
the density of points is proportional to the absolute value 
\q{x)\ while red (green) points mark the regions of posi- 
tive (negative) topological charge density. Note that the 
empty areas correspond to indeed either vanishing or at 
least utterly small 

Several qualitative comments are now in order. It 
is apparent and remarkable that the topological charge 
density is distributed in localized regions which roughly 
look like four dimensional balls. As a matter of fact this 
regions are almost spherically symmetric resembling in- 
stantons from the first sight. We are in hustle to add how- 
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ever that we don't claim that these regions are indeed the 
instantons (see below) . In either case the lumpy structure 
of the topological charge density is present in all configu- 
rations we have and could be taken as firmly established. 
Note that this result is not new, the lumpy distribution of 
the topological charge in "hot" (not cooled, smoothed or 
smeared in any way) vacuum configurations was observed 
already in Refs. \l4. Il5l Il6l Il7| . However, we found the 
lumps completely independently using newly developed 
method which not only confirms their existence but also 
provides new opportunities to investigate the structure of 
the lumps. The second point is that the lumps Ci in q(x) 
distribution are significantly dilute, their total number in 
every time slice could be estimated as < 10, the corre- 
sponding lumps density being Dc ~ 3 • 10~ 3 in lattice 
units which translates to Dc ~ 10 fm~ 4 for (3 = 2.40. At 
the same time it is apparent that the lumps are clearly 
separated from each other and their characteristic size is 
of order p c « (1 -e- 2) a « 0.12 0.24 fm. 

The reminder of this section is organized as follows. In 
section we describe in details our approach aimed 

to investigate the lumpy structure described above. In 
particular, we argue that in order to quantitatively de- 
fine the notion of lumps it is mandatory to introduce the 
physically motivated cutoff A q on \q(x)\ values and to en- 
sure that final results are cutoff independent in reason- 
ably large range of A q . In section II V C 21 we explore the 
limit A q — > which seems to be unphysical for us, but 
in which we find that the topological charge organizes 
into some global structures reminiscent to ones discov- 
ered recently 0. Then in section HV C 31 we present our 
results for the topological charge density correlation func- 
tion and discuss some related issues. Until section llV C 41 
we consider only the data collected on 16 4 , j3 = 2.40 lat- 
tices, while in section lTV C 4l the comparison is made with 
the results obtained at (3 = 2.475. 



1. Lumps in q(x) Distribution 

In order to investigate quantitatively the nature of the 
lumps Ci in the topological charge density we have to de- 
fine first an algorithm to identify the lumps. The natural 
procedure is to apply some cutoff A q to the topologi- 
cal density and to treat the small values \q{x)\ < A q as 
identical zero. It is expected that for reasonably wide 
range of A q the lumps £, would look like an isolated 4- 
dimensional regions of sign-coherent topological charge. 
It is clear that the number of lumps and other their char- 
acteristics depend upon the cutoff A q . We expect that 
the universal physical properties of the lumps in q{x) dis- 
tribution should be independent or depend only mildly 
on A q (in the second case the bias introduced by A g de- 
pendence is to be included into systematic errors). 

Notice the following two extreme cases. Evidently, for 
large A q there would be no lumps at all and hence the 
maximal cutoff coincides by order of magnitude with the 
maximal value of On the other hand for A q = 
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FIG. 13: Dependence of the lumps density Dc upon the cutoff 
A 9 for various values of Vmin, Eq. (I3H . 



one expects to find a large amount of very small lumps 
consisting of just one lattice point. There are at least 
two reasons to justify this assumption. First, this is 
because we constructed the HP 1 a- model starting from 
"hot" thermalized gauge configurations and | q ) fields re- 
flect partially the ultraviolet noise present in the original 
gauge potentials. Therefore, away from the lumps the 
topological density is expected to fluctuate widely around 
zero while being almost vanishing in magnitude. Sec- 
ondly, it is quite evident that it makes no sense to assign 
any physical meaning to utterly small |?(a;)|. The phys- 
ically meaningful minimal Aq value could be estimated 
by order of magnitude as the ratio of typical topological 
charge and the lattice volume, which for our lattices is 
« 10~ 5 . We conclude therefore that the extremely small 
cutoff A q ■ a 4 < 10~ 5 could not have any physical sig- 
nificance and the universal properties of the topological 
density bulk distribution (universality window) should be 



looked for in A„ 



> 10 region. 



Let us consider first the dependence of the total lumps 
density Dc upon the cutoff A g , which is shown by the 
upper curve on Fig. 1131 As expected from the discussion 
above Dc is monotonically decreasing function of A q and 
diverges in the limit A, -> 0. However, the dependence 
Dc(A-q) becomes less trivial if we consider only lumps 
satisfying 



V(d) > V mm >l 



(31) 



where V(Ci) is the total number of sites constituting the 
lump Ci and V m i n is taken to be of order few units. In 
fact, the constraint (|31|l is quite natural since for V m i n = 
1 it cuts only smallest fraction of UV-scale lumps. It is 
remarkable that already from V m in — 1 onward the shape 
of Dc{A q ) dependence is universal, the only thing which 
distinguishes the curves at various V m i n is the scale of 
the density Dc(A q ). Note that the behavior of Dc in the 
limit Aq — > is changing qualitatively once the volume 
cutoff is imposed. As is apparent from Fig. 1141 for any 
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FIG. 14: Dc(kq) dependence in the range 10~ 5 < A q ■ a 4 < 
1.5 • 1(T 3 for various V min , Eq. 




FIG. 15: Distribution of lumps volumes, Eq. (1341 . for A q -a 4 — 
1.5 • 10 -3 calculated on 16 4 , /3 = 2.40 lattices and normal- 
ized by the number of configurations. Solid line is the fit to 
Eq. Eg). 



Vmin > we have 



lim Dc 



V mi „>0 



const < oo 



(32) 



and this is discussed in the next section. The maximum 
of the function Dc{A q ) at V m i n > and its rather mild 
dependence on A q for 



A, > 1(T 3 ■ a' 4 = (290 MeV) 4 



(33) 



suggests that A g universality window starts from this 
value. 

Let us now consider the volumes and the topological 
charges of the lumps £i 

v{d) = Yl 1 > = E ' ( 34 ) 




FIG. 16: Lumps charge, Eq. 13411 . distribution for A g -a 4 = 1.5- 
10 -3 calculated on 16 4 , (3 = 2.40 lattices and normalized by 
the number of configurations. Solid line is the fit to Eq. 1351 . 



identified for given value of A q . The distributions of 
V(£i) andQ(A) for A,-a 4 = 1.5- 10" 3 calculated on 16 4 , 
(3 = 2 AO lattices and normalized by the number of con- 
figurations are shown on Fig. 1151 and Fig . 1 1 61 respectively. 
A few points are worth to be mentioned here. First, we 
see that both volume and charge distributions are well 
described by power laws 



N[V(d)] = c v ■ [V{CA]- av , 
N[Q(d)} = c Q ■ [Q(A)r QQ 



(35) 



for not exceptionally large V(£i) and Q(A)- O n the 
other hand the appearance of extremely large lumps with 
volumes V{Ci) > 100 • a 4 and charges Q{£i) > 1 is 
cutoff dependent (see below). Therefore, let us concen- 
trate on the behavior (|35|l valid for moderate V(Ci) and 
Q(Ci). It turns out that both histograms in the almost 
entire ranges of volumes and charges are well described 
by Eq. (|35|l and the fitted values of power exponents are 
given by 



"Ha, -^=0.0015 = 1-52(8), 
a Q I a„ .„«=o.ooi5 = 1-32(20), 



(36) 



where a conservative error estimates are quoted. Ap- 
parently both power exponents are compatible with each 
other and it is tempting to conclude that for A g • a 4 = 
1.5 • 10~ 3 they are equal to 



ay = ctQ = 3/2. 



(37) 



Clearly, we should check this conjecture at very least by 
investigating the dependence of power exponents upon 
the cutoff A q . The corresponding histograms are pre- 
sented on Fig. El an d Fig. ED where we have collected 
the distributions of V(Ci) and Q(£i) obtained for vari- 
ous A g on 16 4 , (3 = 2 AO configurations. The solid lines on 
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FIG. 17: Lumps volume distribution for various A q given in 
lattice units. Solid line is the fit 135H with ay = 3/2. 



FIG. 19: Lumps limited squared topological charge (Q 2 )i um p 
versus cutoff A 9 . Lines represent {Q 2 }, Eq. 13011 . and the 
corresponding error bounds. 




FIG. 18: Distribution of lumps charges for different A q given 
in lattice units. Solid line is the fit 135H with olq = 3/2. 



the figures represent the one-parameter fits to Eq. Ij35(l 
in which we have fixed ay and cxq according to Eq. (|37|l . 
It is remarkable that the volumes of the lumps are dis- 
tributed well within numerical uncertainties in agreement 
with I35|) at all values of the cutoff. The only deviation 
from l|33|l is seen at large A q , V(£i) and could be easily 
understood. Indeed, for very large A q we will not find 
any lumps at all and with diminishing cutoff we start 
seeing lumps from some finite volume only. This could 
explain the apparent drop in the V(Ci) distribution for 
A q ■ a 4 equal to 6 • 10 -3 and 10~ 2 . Essentially the same 
arguments apply to the Q{C{) histograms which show a 
slight deviation from (|35[l for large A q and small Q(Ci). 
Thus the have found an universal (cutoff independent) 
behavior (|35() , (|37|l of the lumps volume and charge dis- 



tributions for A„ 



> 10- 



In the next section we 



show that for smaller A„ this universality does not hold. 
Therefore we have clear evidences that the lower bound 



of the cutoff universality window is indeed given by 13311 . 

Let us make a few comments concerning the distribu- 
tion l|35|) . (|37|l . The simplest possible picture behind the 
lumpy structure of the topological density is that q(x) is 
non-vanishing only within the four dimensional balls of 
radius p. It is reminiscent to the well known instantonic 
picture of the vacuum and in this case the power law 
(|35|l . (|37|l translates into 



N[p]dp ~ ^ 



dp 



(38) 



which is precisely the instantons distribution elaborated 
in Ref. |33j. While the microscopic view of the lumps as 
four dimensional balls of sign-coherent topological charge 
looks quite natural, it is in fact incompatible with instan- 
tons. Indeed, Eqs. (|35|l . (|37|l imply that the lumps are 
objects with fractional and arbitrary small charge. As far 
as the large lumps with Q{d) sw 1 are concerned we dis- 
cuss them in the next section, but right now it is enough 
to mention that their number is negligible compared to 
the total number of lumps. The conclusion is that the 
majority of the lumps are not the instantons. 

The question of primary importance is whether the to- 
tal topological charge is dominated by the lumps. The 
relevant quantity is the lumps limited topological suscep- 
tibility which is essentially (Q 2 ) measured on the lumps 
only. Note that this quantity naturally bounds the physi- 
cally acceptable cutoff from above. Indeed, for vanishing 
cutoff the lumps limited average (Q 2 )i U m P coincides with 



the full one, (Q 



Hump 



(Q ), while at very large A q we 



would have (Q 2 )i um p = 0. The characteristic cutoff at 
which the topological susceptibility drops down is defi- 
nitely the upper bound on physical A q and a priori could 
be on the either side of the value quoted in (|33p. The 
lumps would be clearly unphysical if A q range happens 
to be empty. 
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We have measured the lumps limited squared topolog- 
ical charge (Q 2 }i ump on 16 4 (3 = 2.40 configurations, the 
result is presented on Fig. EH Note that we also checked 
the dependence of (Q 2 )i um p on the volume cutoff 131|) 
which turns out to be trivial as is clear from the figure. 
Indeed, the data points with volume cut V m in = 4 im- 
posed are falling on the top of the data with no volume 
cut at all. It is apparent that the topological susceptibil- 
ity is indeed lumps dominated provided that 



A < 2 



l(T 3 a- 4 



(350 MeV) 4 



(39) 



Note that the universality window (|33|l . I|39|) appears to 
be rather small. However, this does not appear com- 
pletely unexpected. Indeed, our original intent was to 
find a plateau in the cutoff dependence of various quan- 
tities characterizing the lumps, at which their first A g 
derivative vanishes. In particular, the universality win- 
dow (|33|l . (|39|l could be estimated rather precisely by 
noting that the maximum of the lumps density Djc(A q ) 
at any Vmin ^ 1 is located roughly at 



A* 



(lH-1.5) 10" 



(314(16) MeV) 4 , (40) 



which almost coincides with the universality window 
range. 

Finally, let us note that the question of dislocations, 
which a priori could be crucial for our investigation, is 
likely to be irrelevant in fact. For general arguments, 
which we could just verbosely repeat here and which show 
that dislocations are not essential in our case, the reader 
is referred to Ref. [l7|. We could only add that the dis- 
location filtering (as well as partial filtering of the UV 
noise) is built into our approach. Moreover, we don't 
see any sign of dislocations in the topological susceptibil- 
ity which seems to be perfectly scaling and is in a good 
agreement with its accepted value. Following Ref. 0] 
we stress that the notion of dislocation is indispensable 
from both the gauge action and the topological charge 
operator. Therefore due to the specific properties of our 
topological charge construction we expect that the issue 
of dislocations inherent to Wilson gauge action is not ap- 
plicable in our case. 



2. Exploring the Limit A q — > 

In this section we explore the limit of vanishing cutoff 
imposed on the topological charge density. Although it is 
not clear for us what is the physical relevance of the limit 
A q — > it nevertheless seems to be interesting to consider 
in connection with recent observation of low-dimensional 
global structures in the topological charge density dis- 
tribution [l|| . Essentially, these global structures are 
defined as topological charge sign-coherent regions and 
thus are similar to the lumps we're investigating. 

In fact, the non-trivial features of the A g — > limit 
could be seen already on Fig. ^] from which it follows 
that once the mildest V m in = 1 volume cut is imposed 
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FIG. 20: Lumps volume distribution for fixed cutoff A q ■ a — 
2-l(T 5 . 
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FIG. 21: Volumes of percolating lumps at various A q . Note 
that percolating lumps disappear at A q ■ a 4 w 7 • 10 -4 . 



the density of lumps rapidly diminishes with decreasing 
A g provided that A g -a 4 < 10^ 3 . The dependence Dc(A q ) 
in the range (see Fig. I14fl 



A q e [2 • 10" 



io- 



(41) 



= [(110 MeV) 4 ; (290 MeV) 4 



makes it apparent that Dc(0) is very small if not non- 
vanishing. However, the total number of lumps (no vol- 
ume cut imposed) is divergent, hence the dominant con- 
tribution at small A q comes from lumps consisting of just 
one point. Then the question is whether we have any 
other contributions to Dc or everything is exhausted by 
one-point lumps. 

To investigate this issue we considered the lumps vol- 
ume distribution at A q ■ a 4 — 2 • 10 -5 and it turns out (see 
Fig. I20|) that it is qualitatively distinct from the case of 
moderate cutoff. First, we clearly see that the histogram 
separates into two disconnected pieces corresponding to 
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FIG. 22: Distribution of the small volume lumps at various 
A 9 . Note that it still follows the power law 135H . 

small and extremely large (percolating) lumps. As far 
as the percolating lumps are concerned we didn't in- 
vestigated their properties in details for reasons to be 
explained shortly. We only mention that there are al- 
most always only two percolating lumps on our configu- 
rations carrying an extremely large and opposite topolog- 
ical charge. Qualitatively our results are in accord with 
that of Refs. jig, where the analogous percolating global 
structures were found. If we now turn to the consider- 
ation of small lumps distribution, it also follows sharply 
the law (|35[) . However, the power exponent is drastically 
changing 

«v| Ag . a 4 =2 . 10 -5 =2.9(1). (42) 

Next we compare the lumps volume distributions at 
various A q from the range (|41|l . The first observation 
is that the volume of the percolating lumps strongly de- 
pends upon the cutoff imposed as is shown on Fig. 1211 
Moreover, the percolating lumps are present only for 
A g • a 4 < 7 TCP 4 and disappear at larger values of the cut- 
off. At the same time the power exponent ay is also not 
universal as is clear from Fig. 1221 where we have plotted 
the distribution of lumps with small volumes. We have 
also measured the function ay(A q ) and it is presented 
on Fig. 1231 The power exponent ay becomes cutoff inde- 
pendent only for A q ■ a 4 > 1CP 3 thus justifying the choice 
(1331) . I|39l) made earlier for the physical range of the cutoff 
values. 

To summarize, it seems crucial to ensure the univer- 
sality of the lumps distribution (and thus the very def- 
inition of the lumps) with respect to the a priori arbi- 
trary choice of the cutoff on the topological charge den- 
sity. Once the universality is ensured the lumps distri- 
bution follows sharply the cutoff independent power law 
(|35|l with rather remarkable value of the power exponents 
ay « ojq w 3/2. To the contrary at small A q we found 
very peculiar lumps structure. Namely the number of 
small lumps obeys the same power law Il35f but with 



FIG. 23: The power exponent ay, Eq. 1351 . as function of A 9 
in the range I4H . Line represents the conjectured universal 
value JSZJ' 

much large power exponent which seems to be divergent 
in the limit A q — > 0. At the same time a few (typically 
two) percolating lumps appear which are reminiscent to 
the global topological structures discovered recently [18| . 
However, it is unclear for us what is the physical signif- 
icance of the results obtained at very small A q since the 
lumps properties are strongly cutoff dependent in this 
limit. 



3. Topological Density Correlation Function 

In this section we present our results for the topologi- 
cal charge density correlation function (q(0)q(x)) = (qq) 
obtained on our 16 4 , (3 = 2.40 configurations. However, 
let us clarify first our view on some theoretical points 
concerning this correlator which were much elaborated 
in the recent studies [30(. 

It has long been known (see Ref. |34j and references 
therein) that the correlation function (qq) can not be pos- 
itive at non-zero physical distances. Confronted with the 
positivity of the topological susceptibility x — (Q 2 )/V = 
J(qq) this implies a peculiar structure of (qq) correlation 
function, which is to be negative for x 7^ 0, while con- 
taining non-integrable singularity (contact term) at the 
origin such that the susceptibility remains positive (in 
the context of CP™ models this issue is also discussed 
in Ref. H^)- In fact the negativity of (qq) correlator 
motivated the recent discovery of the global topological 
structures since it implies that the topological charge 
could not be concentrated in four-dimensional regions of 
finite physical size. However, we do think that the con- 
tact term in (qq) correlator is the clear sign of the per- 
turbation theory to which the susceptibility is believed 
to be unrelated. 

To sort out the problem consider the instanton liquid 
model where the identification of perturbative and non- 
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perturbative contributions is intuitively clean. At least 
at small non-zero distances we have (qq) per t ~ — 1/ 1 a; | 8 , 
while {qq) non-pert is positive at small x and rapidly drops 
down at distances comparable with the characteristic in- 
stanton size (it could even become negative at larger x if 
the instantons interaction is significant). If we suppose 
that the perturbative and non-perturbative parts do not 
mix then the negativity requirement could be fulfilled at 
least in principle. The outcome is that the perturba- 
tion theory could overwhelm the non-perturbative part 
and reproduce qualitatively the behavior of (qq) outlined 
above still giving vanishing contribution to the suscepti- 
bility. Then the actual question is how do we understand 
the topological charge density operator. The negativity 
requirement could hold for perturbation theory sensitive 
(thus presumably local) definitions of q(x) and does not 
apply to non-perturbative ones. Thus the problem we're 
discussing is essentially the same long standing problem 
of perturbative and non-perturbative physics separation, 
which seems to have no unambiguous solution (see, e.g., 
Ref. (361] for recent discussion). 

A particular illustration of the above reasoning is pro- 
vided by the overlap-based construction of the topolog- 
ical charge (lC||. It is remarkable that once only a few 
lowest modes of the Dirac operator are considered the 
topological charge density distribution has lumpy struc- 
ture 0, 0, 0| while the (qq) correlation function stays 
non-negative at all distances. However, the inclusion of 
higher Dirac modes into the q(x) definition results in 
the appearance of some global "sheets" percolating 
through all the lattice volume and carrying extremely 
large topological charge. Simultaneously with increas- 
ing number of modes the (qq) correlation function be- 
comes negative and develops the positive singularity at 
the origin (30f. The interpretation might be that the 
more modes are included the more local the definition of 
topological density is and hence q(x) becomes more sen- 
sitive to the perturbation theory. Essentially for this rea- 
son we considered the effective topological density (|25[) 
in section ITlIB 21 where the overlap-based definition was 
compared with our construction. 

It is remarkable that we also see qualitatively the same 
behavior with our definition of the topological charge 
density. The border between perturbative and non- 
perturbative regions is provided in our case by the uni- 
versality window with respect to A q . We have seen that 
the lumpy structure gradually changes with diminishing 
A q and becomes consisting of a few percolating lumps 
plus the divergent number of lumps having UV-scale vol- 
ume. However, we stress that not only the window (1331) . 
(|39[) provides the separation of perturbative and non- 
perturbative physics. Some degree of non-locality is built 
into our approach from very beginning and in the spirit of 
the above discussion one could say that our definition of 
the topological density is non-perturbative by construc- 
tion to which the negativity requirement does not apply. 
Thus we expect that the (qq) correlation function in our 
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FIG. 24: The topological charge density correlation function 
measured on 16 4 , /3 = 2.40 lattices for A q -a 4 = and 1.5-10~ 3 . 
Lines represent the best fits 1431 . 



case should behave roughly exponentially 
(qq) ~ exp{-\x\/p c } 



(43) 



with correlation length p£ corresponding to the charac- 
teristic lump size. 

The result of our measurements of the (qq) correlator 
on 16 4 , j3 — 2.40 configurations is presented on Fig.^J It 
should be noted that within the numerical errors (qq) is 
nowhere becoming negative and is compatible with zero 
from |x| = 6 • a onward. Therefore we plotted the (qq) 
correlator only in the range \x\/a e [0; 5] where data 
closely follows Eq. (|43|l . In turn, the fitted value of the 
characteristic lump size is given by 

Pc\ Aq=0 = 0.50(2) -a = 0.060(3) fm, (44) 

where we have explicitly indicated that the data was 
taken without any cutoff imposed. As far as the depen- 
dence upon the cuts (J2HJ), and (|3*T)l is concerned, it 
turns out that the (qq) correlator is independent on V m i n 
at non-zero separations provided that V m i n is of order 
unity. Indeed, the dominant contribution to (qq) comes 
from the extended regions of sign-coherent topological 
charge while the contribution of small lumps V(Ci) ~ 1 
averages to zero for x ^ 0. Note that the apparent devi- 
ation of (qq) from l|43l) at \x\/a = 0, 1 is due to the abun- 
dance of small lumps and practically disappears once the 
cut Vmin ~ 1 is applied. The dependence upon A q could 
also be qualitatively established since for increasing cut- 
off the size of the lump contributing significantly to (qq) 
diminishes. In fact, pc depends linearly on A q and for 
A q ■ a 4 = 1.5 • 10~ 3 it is given by 

Aclva*=i.5.io-> = 0.41(2). a = 0.049(2) fin. (45) 

To summarize, the topological charge correlation func- 
tion (qq) respects the above investigated lumpy structure 
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of the topological charge bulk distribution. Namely, it 
falls exponentially with distance, the corresponding char- 
acteristic size of the lumps is given by B44JI and decreases 
if the non-zero cutoff A q is applied, Eq. H45fl . The quoted 
pc values are in qualitative agreement with lumps vol- 
ume distribution we investigated previously, namely, it 
is rather small in lattice units respecting the ultraviolet 
divergence present in (|57l) . 

4- Scaling Check 

In this section we confront the data obtained at (3 — 
2.40 and (3 = 2.475 in order to check the scaling proper- 
ties of the quantities introduced earlier. First, consider 
the physical range of the cutoff A,j which could be char- 
acterized by A*, Eq. (|4*U)l . Performing the analogous cal- 
culations for (3 — 2.475 we found that 

A* = (4 8) • 1CT 4 • a~ 4 = (340(28) MeV) 4 , (46) 

which is consistent with its value at (3 = 2.40. Thus 
the cutoff introduced to define the lumps appears to be 
physical quantity with continuum value around (300 ~ 
350 MeV) 4 . Note that the range of A, where d Aq D c 
approximately vanishes is notably wider at [3 = 2.475 
then is was at [3 = 2.40. 

As far as the power exponents ay and «q are con- 
cerned, at (3 — 2.475 they turn out to be also compatible 
with each other. In particular, ay values at the bound- 
aries of the range (|46[1 are given by 

«vH A9 . a 4 =4 . 10 -4 = 1-62(14), (47) 
"v| Ag . a 4 =8 . 10 -4 = 1.56(12) 

and are consistent with the conjecture l|37|) . On the other 
hand, the topological susceptibility is lumps saturated for 
A g -a 4 < 10~ 3 which is even a bit beyond the upper bound 
quoted in i|46|) . Therefore the estimate of the universality 
window (|46[) is in fact rather conservative. 

Finally, we considered the topological density correla- 
tion function at (3 — 2.475. It turns out that it is again 
non-negative within the numerical errors and falls off ex- 
ponentially with distance. The corresponding character- 
istic lump size at j3 = 2.475 and A,j = is 

Pc\ Aq=a = 0.57(2) -a = 0.052(2) fm, (48) 

which disagrees with (|44|l and shows instead that 
Ac(A g = 0) is almost constant in lattice units. How- 
ever, this should not be surprising since A, = is def- 
initely outside the physically acceptable region. On the 
other hand, the characteristic lump size measured at 
A q ■ a 4 — 6 • 10~ 4 turns out to be 

AclA..a'=6.io-' = 0.50(1)- a = 0.046(1) fm (49) 

and is compatible with l|45[) . At present it is too early to 
speculate about the observed scaling in the topological 
charge density correlation function, obviously more data 
is needed to quantify the issue. 



V. CONCLUSIONS 

In this paper we considered the alternative definition of 
the topological charge density in pure SU(2) lattice gauge 
model. Generically the idea is to exploit the well known 
close connection between HP" cr-models and SU(2) Yang- 
Mills theory, which was proved to be very useful in the 
investigations of Yang-Mills fields topology. The usual 
approach in the past was to consider the cr-model instan- 
tons and to induce the corresponding SU(2) gauge poten- 
tials which realize the topologically non-trivial configu- 
rations in the gauge theory. We tried instead to read the 
above relation in the opposite way, namely, to find the 
unique HP" cr-model fields which are closest to the given 
SU(2) gauge background. In the continuum limit and for 
the smooth fields we are confident that our construction 
is identical to any other way of the topological charge 
calculation. However, considered within the regularized 
quantum field theory the HP" cr-model approach obtains 
a separate significance and provides an alternative defi- 
nition of the topological charge on the lattice. Moreover, 
the cr-model construction naturally leads to the unam- 
biguous definition of the topological charge density with 
clean geometrical meaning. 

The actual realization of the above program brought 
out a wealth of technical issues which, we believe, were 
adequately addressed in the paper. Moreover, the algo- 
rithm was tested in various circumstances and, in par- 
ticular, it was compared with overlap-based topological 
density definition. In the latter case we found a strong 
evidences that both approaches are consistent with each 
other. 

As far as the topology of SU(2) gauge fields are con- 
cerned the results we obtained are at least in qualitative 
agreement with the literature. In particular, the topo- 
logical susceptibility measured by our method seems to 
be perfectly scaling and coincides with its conventional 
value. Note that at present there is no commonly ac- 
cepted picture of the topological charge bulk distribution. 
Here our approach provides the unique investigation tool 
and our main results could be summarized as follows. 

We confirm the lumpy structure of the topological 
charge bulk distribution discovered in the previous stud- 
ies and assert that these lumps are not compatible with 
instantons. Moreover, for carefully chosen parameters 
entering the lumps definition we found that their volume 
distribution seems to obey rather remarkable power low 

N(Vc) ~ Vr 3 ^ 2 , while their topological charge is likely 
to be proportional to the lumps volume with distribu- 
tion being essentially the same. At the same time the 
lumps do saturate the topological susceptibility and in 
the first approximation the topological charge is indeed 
concentrated in the lumps only. 

Note that the lumps volume distribution naively im- 
plies that the characteristic size of the lumps is given by 
the UV scale. On the other hand, it could be probed 
by the topological density correlation function (qq) pro- 
vided that we could separate the perturbative contribu- 
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tion from it. We argued that our topological density 
construction is inherently non-perturbative and its cor- 
relation function generically does not obey the negativ- 
ity requirement. In turn the data indicate strongly that 
the (qq) correlator falls off exponentially, the correla- 
tion length being indeed of order half the lattice spac- 
ing. However, it increases with diminishing spacing so 
that the possible scaling is not excluded. In cither case 
the characteristic lumps size extracted from (qq) correla- 
tion function is definitely smaller then 0.1 fm « 2 GeV -1 
and is worth to be compared with the approximate loca- 
tion of A„ universality window (300 350 MeV) 4 . The 
scale of few hundred MeV is inherent the gauge fields 
topology, for instance, it is built into phenomenological 
instantonic models. In our approach this scale is the 
lower bound on the characteristic topological density in 
the lumps and separates the lumps from the UV noise 
present in q(x). Thus at zero approximation the lumps 
could be viewed as highly localized pc < 2 GeV -1 and 
'hot' qc > [300 MeV] 4 bumps in the topological density, 



where pc and qc are the characteristic lumps size and 
density correspondingly. Therefore they could strongly 
influence the microscopic structure of low-lying Dirac 
eigenmodes. In particular, the lumps smallness could 
explain the unusual localization properties of low Dirac 
eigenmodes discovered recently 0, l37l l38| , while their 
characteristic topological density could reveal itself in the 
Dirac modes mobility edge [38| . 
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